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The singular integral equation 



1 
1 2 



U(x) cfo=A, 0<(7<1, 



x— a x-\-<jj 
is converted into another singular integral equation 

2ji t— x 2 v Ji s/t(t—x) ^x 

where 



<j>(x)- 



L'x' 



* ' IxJ 



■\/x 

The former was derived from a boundary value problem in wave guide theory, and the solution is known. By 
solving a related Hilbert problem for two unknown functions, the solution of the latter integral equation is 
obtained. Then, the expression 



*(*)-i*(?) 



is compared to the former known expression obtained for \p{x), and the new expression is found to be more 
tractable. 

Key words: Analytic function; Cauchy principle value; Hilbert problem for two unknowns; singular integral 
equation; wave-guide boundary problem. 

1. Introduction 

This paper illustrates an alternate procedure for solving in closed form an integral equation 
derived from a problem concerning a 3-dimensional wave-guide structure and discussed by 
G. T. Bierman [l] 1 namely 

/ 1 \ 

t(x)dx=h, 0<cr<l. (1) 

The more general class of equations, of which (1) is a special case, in also treated by H. F. 
Bueckner [2] using numerical approximation methods. 

One purpose for the alternate procedure is to arrive at the solution is a more simplified 
expression than previously obtained. Another purpose is also intended. The method we employ 
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involves solving a Hilbert problem for two unknowns. Although the method for solving the 
general Hilbert problem for one unknown is known, [3], [5], [7], Vekua [10] states that the 
solution to the general problem for several unknowns cannot be expressed in closed form. 
The other purpose is to illustrate a special technique which can be used to solve certain non- 
homogeneous Hilbert problems for several unknowns. 

2. Physical Origin of (1) 

Equation (1) has been shown by Bierman [1] to arise from a problem involving a wave- 
guide structure in which a perfectly conducting surface extends from — oo to oo in the direc- 
tion of the z axis, making the boundary problem 2-dimensional in nature. Although some 
numerical results for the physical problem are reported in [6] and [9], the objective in [1] is to 
obtain a closed form solution. Here, we summarize the results beginning with the following 
diagram taken from [1] of the boundary and boundary conditions in the x-y plane. 



27T <*- 
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Figure 1. Widening, straight, 2 -dimensional wave guide 



E e ikx 



By Fourier analysis, the following equation for 0<y<qir is obtained. 

0-) 

\£/ , g sin (g) 



^-'^(a-I^Hf *' tt) 



sin i 



cos 



(!)-(!) 



cos (£)— cos (y) 



df, 



(2) 



where E(x, y) is the ^/-component of the electric field of the wave guide, and E(y)=E(0, y) 
represents the field in the aperture. Solving for E'(y) and hence E(y) then determines the 
reflection coefficient p by 



Once p is computed, the transmission coefficient is given by 

T=q(l-p). 

Although the general solution to (2) for arbitrary g, 0<g<l, is unknown, Bierman arrives 
at a solution for the specific case 

1 
2=2 • 



(2a) 



(2b) 
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This is done by letting 

z=cos (y), a=cos ({), +(x)=E'(cos- 1 (x)), h=^ik(l-p)E (^-^^ (2c) 

in which case (2) for 



becomes 



r 



<x — a X-\-a 



*=2 



— h(x)dx=h, 0<cr<l. (1) 



Then (1) is solved for \p(x) and E'(y) obtained from (2c). 

Actually, a general class of equations of which (1) is a specific case is treated in [1]. In 
the treatment, the general equation is converted into an equivalent integral equation to which 
the Wiener-Hopf technique is applied, i.e., the technique of Fourier transforms, functional 
equations, and analytic continuation procedures. The general solution obtained for (1) is 

^)=^)+g|J o 7>„(f)i 1 (|)i^, (3) 



where 



and 



1 — 2/3 

£„(*)= T ^==[(l_VT=^) 2/3 +(l+VT=^) 2/3 ], 0<*<1, (3a) 



i^Tk](x)= 



|^r 0<i<1 ' 



-4 J I cosh (^- £j+3 ( i±rt) sinh (3 y\ 7t2Z_lp d t> w=x-jx 2 —l, x>l. 



Expression k (x) of (3a) is a nontrivial eigenf unction of (2), and c is a constant. 

A disadvantage of the above result is the complicated expression for the particular part 
of the solution to (1). Limitation due to the difficulty in evaluating ki(x) and the particular 
part of the solution (3) is admitted in [1] by stating that it does not seem likely they can be 
expressed in terms of elementary functions. However, we shall show that the solution can be 
expressed in terms of elementary functions. 

3. An Alternate Approach to Solving (1) 

In this section, we shall convert (1) into a singular integral equation in which only the 
Cauchy kernel term j— appears under the integral sign. This new integral equation will 
subsequently be converted into an equivalent Hilbert problem for two unknown functions [10]. 

By the transformations, 

1 1 
x=->a=-> 
r s 



eq (1) becomes 



r~Kr) 1 s r»Kr) 1 

Ji r r—s 2ji r r+s 
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where we have used a bar through the first integral sign to denote the Cauchy-principle 
value of the integral. Then the substitutions s=^j x , r=^jt, transform (4) into 

r _P <t>(t) dt V» f <j>(t) dt__ h , ,_ , 



where 



Now if the term 



v^ _r *w ii 

2 Ji (y/t—y/x) 2-Jt 



is added and subtracted on the left-hand side of (5a) , then 

Also, if the term 

0(0 ^ 



^r 



(5b) 



(6) 



(Vi+V*) 2 V< 

is added and subtracted on the left-hand side of (5a), then 
Finally, by adding eqs (6) and (7) , 

irm dt+ irr-_ML. dt== ^ n< ^ (8) 

2 Ji £ — x 2 Ji ^t^t — x) ^x 

where 

* (a;) -^r' (9) 

Equation (8) is a special case of the singular integral equation 

a(x)-£ m t^dt + bix)^ r^ dt=c(x), l<z<> (10) 

Ji * — * Ji -yjt(t-x) 

where the bar through each integral sign denotes the Cauchy principle value of the integral. 
This class of equations is considered in [8] where the approach is to consider a Hilbert problem 
derived from (10). Another case of (10) has been shown by Latta [4] to arise from a convolution 
integral equation with finite limits. In the next section we begin discussing the Hilbert problem 
related to (8) . 

4. Conversion of (8) to a Hilbert Problem for Two Unknowns 

To solve (8) , one can try to proceed with a method similar to one of the standard approaches 
used to solve certain singular integral equations with the Cauchy kernel term j— - [3], [5], [7]. 
The method is based on the theory of Cauchy integrals [5], [7]. 
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If a solution is known to exist, the theory of Cauchy integrals suggests its derivation (for 
reasons which will become clear in the ensuing analysis) by considering the functions 

If <f>(t) is of a suitable class of functions, 2 then F(z) and G{z) have the following properties: 

I. F(z) and 6?(V) are analytic in the finite complex 2-plane excluding the real semi- 
infinite interval [1, 00). 
II. F(z) and 6{z) have at worst growth rate at z=\ less than that of a simple pole. 
III. There exists a y>0 such that 

F(z)=0(\z\-y) 
and 



G(z) 



= o(|s| V 7+ 2)) as | |_>go, ze[l, 00). 



Other basic properties of F(z) and 6? (2) come from relations between their boundary values. 
Let F+(x), G+(x) and F_(a;), #_(;£) denote the limits of F(z), G(z) as z— >x, 1< #< 00 , from above 
and below the real axis respectively. From the theory of Cauchy integrals [5], [7], these limits 
exist and are given by 3 

which leads to the Plemelj formulae 

F + (z)+^_(z)=-i.-f"^A (13a) 

TTl Jl t — X 

F + (x)-F_(x)=<j>(x), (13b) 

G+{x)+G_(x)=X-f° ,- m dt, (13c) 

G+(x)-G-(x)=^> (13d) 

l<a;< co. 

Now, the integral equation (8) can be expressed as an equivalent problem in terms of 
F(z) and G(z) denned in (11). Note that (13a) and (13c) reveal that (8) is equivalent to 

^Ti(F+(x)+F.(x))+^ 1 ri(6 + (x)+G-(x))=^>l<x<«> l (14a) 



2 It is sufficient that — V^(V^) be of class C(A) described by Levinson [5] where C=(Q, 1). Then, by 
results in [5], the function 



m,(z)=-L f 1 : 

2iriJo 



t-z 



-dt 



is analytic in the complex z-plane cut along [0,1], and has at worst growth rate at z = and z = l less than 
that of a simple pole. Thus, by replacing 

z by - and t by - in y(z), 

z t 

conditions I, II, and III follow for the function 






3 Again, it is sufficient that — vW(Vz) be in class C(A) described by Levinson [5] where C=(0, 1). 
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while (13b) and (13d) reveal that 



F+(z)-F-(x) = y/z(0+(z)-G-(x)), l<Z<oo. 



(14b) 



Equations (14a) and (14b) along with conditions I, II, and III constitute a Hilbert problem 
for two unknown functions F(z) and G(z) [10]. Solution of the problem determines <f>(x) to (8) 
by either (13b) or (13d). In matrix form, (14a) and (14b) can be written as 



c 







ri -3 r -i 

2 ~2^ X 
— 1 -1 






r i 


F + (xy 
+ {x). 


— 


_0-(x)_ 


1 f~2h\ 
~iri \ Jx / 


2 

1 






[_2^x 2 J 






l_ 2 Vz 



!<*<«>, 



which is the matrix analogue of the Hilbert problem on an arc for one unknown function [5], [7]. 
As mentioned in the introduction, the method for solving the general Hilbert problem for 
one unknown is known. However, as Vekua [10] states, the solution to the general Hilbert 
problem for several unknowns cannot be expressed in closed form. One of two special methods 
developed in [8] will be used to solve the Hilbert problem (14). 

5. Solution of the Hilbert Problem (14) 

Relations (14a) and (14b) do not uncouple into equations involving just boundary values 
of F or just boundary values of G. To solve this system, a special procedure which generates 
particular solutions to certain nonhomogeneous problems from known homogeneous solutions 
will be used. 

We first observe that the eigenf unction (3a) to (1) admits an eigenf unction to the cor- 
responding homogeneous form of (8). Specifically, by the relation (9) and the expression (3a), 
the function 



■yJX \X—l 

must be an eigenfunction and the only nontrivial one satisfying 



2Ji t—x 2 Ji Jt(t — X ) 



2 Ji t—x 2 Ji jt(t-x) 
Likewise, by the Plemelj formulae (13) and the integral definitions (11), the function 

V,(s) 

Wis J I I* 



(15) 



(16) 



G Q (z) 



•j_ rm dt 

2wiji t—z 



2xi Ji 



■dt 



, Zi[\, oo), 



(17) 



(18) 



2t%Ji ylt{t-z) 
must be an eigenfunction to 

±vi(F 0+ (x)+F _(x))+^iri(O 0+ (x)+G _(x))=0, 

(F 0+ (x)-F _(x))=^ c (G 0+ (x)-O _(x)),l<x<^. 
To generate candidates for a particular solution to (14), we shall also need the eigenvector 

~Fo(z) 
A(z). 
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satisfying the Hilbert problem "adjoint" to (18), i.e., 

| wi (F 0+ (x)+F _(x))+\ V*« (G 0+ (x)+G _ (x))=0, 

F 0+ (x)-F _(x) = ^x(G 0+ (x)-G _(x)),l<x<^. 



(19) 



The Hilbert problem (19), is of course, derived from the singular integral equation adjoint 
to (16), i.e., 



3 C'Ut) j 1 r _f» Ut) iSrSa, 

2Ji t—x 2 Ji yjttf—x) 

by applying the Plemelj formulae for Cauchy type integrals to 



(20) 



F Q (z) 



G (z) 



'±r~^dt 

2ttiJi t — z 

J_ C"jM_ dt 

2ml J i ijt(t — z) 



(21) 



By the same procedure in which (8) was derived from (1), (20) can be derived from the 
integral equation 



(\^+V^-h^)dx=0, 0<er<l, 
Jo \x—a x-\-a/ 



where 



(t> (x) = - 



(7x). 



(22) 



(22a) 



■yfx 

But (22) can be solved by the method in [1] which shows that there is one nontrivial eigen- 
f unction given by 



-1/3 



Hence, 



^ )= v!^ [(1 - VIz:?)1/3+(1+VI3p)1/3] - 



Ht) i 

ux)=-4M=^L=iUx-jx^iy< 3 +Ux+^x^iy% 



V# -yjx — l 



(23) 



(24) 



and therefore, F (z) and G (z) in (21) are known. 

Now, in order to explain a method in [8] which will construct a candidate for a particular 
solution to (14), we are going to illustrate the idea it involves on a simple Hilbert problem for 
one unknown. This will require less mechanics than if it is explained directly on (14). 

Consider the function 

%)=—=> (*1) 

V2 — 1 

which is analytic in the z-plane cut from z=l to z= «> along the real axis. With the choice of 
the branch 

-yJ^l = ^re id/2 } O<0<2tt, 
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h ± (x) = -- T ==,l<x<co (* 2 ) 

±-\Jx— 1 

where h+(x) and h-(x) are limits of h(z) as 2— >#, l<Cz< oo ; from above and below the cut re- 
spectfully. Hence h(z) is a solution to the homogeneous Hilbert problem 

f + (x)+Mx)=0,l<x<co, (*3) 

where j(z) is required to be analytic in the 2-plane cut from 2=1 to z = oo along the real axis, 
to tend to zero algebraically as |s| — > oo y and to have at worst growth rate at 2=1 less than 
that of a simple pole. Also, b h(z) for any constant b is a solution to (*3). 

Whole families of particular solutions to certain nonhomogeneous Hilbert problems can 
be generated from h(z) in the following way. Let z be a point in the 2-plane not on the real 
axis from z= 1 to z= oo . From 

h(z) 



(z-z ) n 
subtract the principle part of the Laurent expansion of 

Hz) , 

(z-Zo) n 
about z=z 0l and define f n (z) as 

/ rr w Hz) h(z Q ) ^ A<*>(g ) 1___ ; 

/,lW (g-g ) w (s-3o) n tt A! (z-g ) n -* l ; 

where h {k) (z ) is &th derivative of h evaluated at z . For n=l, 2, . . ., a countably infinite 
family of functions is constructed. Each member is analytic in the g-plane cut from z=\ to 
z= co along the real axis. 

To see how the functions (*4) form solutions to certain nonhomogeneous Hilbert problems, 
let 

k=\ 

represent a constant coefficient linear combination of functions defined by (*4). If this linear 
combination converges uniformly to a function g(z) at all finite z in the cut plane and to g+(x), 
g~(x) on each respective side of the cut, then g(z) will be a particular solution to the Hilbert 
problem 

Mx)+Mx)=±h k (x), l<z<oo, (*5) 

k = l 

where 

K*)=M*)-f£z£F k=1 ' 2 ' • • ' ( * 6) 



Conversely, if given the Hilbert problem 

m 

/ + W+/-W=St^T^ Kx<oo, (*7) 



(x—zof 



where b k , k=l, 2, . . ., m, are constants, there is a constant coefficient linear combination of 
functions from (*4), i.e., 

m 

X)X*/*(2), (*8) 

fc=l 

which satisfies (*7). The constants X*, k=l } 2, . . ., m, can be determined by substituting the 
boundary values of (*8) into (*7) and equating coefficients of like powers of 

1 



X Zq 
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The generating function scheme illustrated above can be applied to the Hilbert problem (14) 
to find a particular solution. We have as known homogeneous solutions the eigenvector (17) to 
the problem (18) and the eigenvector (21) to the problem (19). There are several ways to gen- 
erate families of particular solutions to certain nonhomogeneous problems with left-hand side 
as in (14). 

By noting the discussion of the scalar Hilbert problem (*7), it seems as though the fol- 
lowing construction of functions will yield a particular solution to (14) : 



P„(2)1 




r^o(o)-| 


[F (z)-\ 




r^o(o)- 




pV(on 


Z 




z 


z 2 




z 2 




z 


Go(z) 




GWO) 


G (z) 




Go(0) 




G7(0) 


z _ 




_ z m 


L z 2 J 




L 2 2 J 




z 



where F (z) and G (z) are defined in (17). 

Another family of particular solutions can be formed as follows: 



F,(z) 
L0.(s)J 



F (z) F o (0) 



F n+1 (zy 

l.G„+i(z). 



■G n (z) 

F n (z) F n (0) 



= 1,2, 



However, direct substitution will verify that no constant multiple of the first vector in either 
of the above families will satisfy the boundary relation (14b). 

An appropriate candidate for a particular solution to (14) actually comes from the follow- 
ing generation of functions : 



F 1 (z 

l_0l(3)_ 



G (z) 

Fo(z) ft(0) 



F n +i(z) 

\_G„+i{z)_ 



'G n {z) 

F n (z) F n (0) 



n=l, 2, 



where F (z) and G (z) are defined in (21). We will demonstrate that a particular solution to 
the Hilbert problem (14) is a specific constant multiple of 



F,{z) 

l_0i(*OJ 



0o(8) 

Fo(z) F (0) 



(25) 



The constant F (0) is described by 



"• <o) =^X 



t 



dt 



where </> (0 is defined by (24). This definite integral is elementary and not difficult to evaluate 
if one begins with the substitution £ = cosh 2 (0). It turns out that 



F (0) = 



V3i 



(26) 



We shall now show that some constant multiple of (25) is a particular solution of (14). 
Observe that by the way F x (z) and G x (z) were constructed, a constant multiple of (25) satisfies 
conditions I, II, and III for a solution to (14). So, all that is needed is to make it fit the boundary 

relations (14a) and (14b). In terms of F (z) and G (z), the Plemelj formulae (13) for F x (z) 
and Gi(z) in (25) are 

F 1+ (x)-F 1 _(x)=^-> 

sjX 



F + (x) +F 1 _(x)=G 0+ (x)+Go_(x), 
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G 1+ (x)-G 1 _(x)=^ (x), 



(27) 



G 1+ (z)+6i_(x)- 



F 0+ (x)+F _(x) 2F (0) 



!<£<!«>. 



First, note that any constant multiple, X, of (25) satisfies the second relation (14b), i.e., 



\(F 1+ (x)-Fi_(x))=\MOi + (x)-G 1 _(x)), !<*<«>. 



(28a) 



In order for 



L<?i(s)J 

to satisfy relation (14a), X must be chosen so that 

^Ti\(F l ,(x)+F 1 _(x))+^TiUGiM+G,_(x))- ' 



l<Cz<°° . 



(28b) 



By substituting the boundary relations for ^1(2) and G x {z) in (27), (28b) becomes 

, Tl . /« ~ x . 3 iti t ~ ~ \1 , 3 r ., /-2#o(0)\ -2A ,„. 

X |_2 « (<?o + (^) + c?o-(x))+2^ (Fo t (^)+F -(a:))J+2 Va; ' rlX ( 3 j = ^"' (29) 

By the known relation (19), (29) reduces to 



|V^x( 



2^ (0) N 



-2A 



V5 



(30) 



Finally, by substituting (26) into (30) and solving for X, 

irV3 
What has just been shown is that the vector 

"^i(2) _ 



h 
xV3 



is a particular solution to the Hilbert problem (14). The general solution to (14) is an arbitrary 
constant multiple of the homogeneous solution (17) plus a particular solution. Therefore, the 
general solution to (14) is given by 



-F(z)i r^>(z)- 

lG{z)_~ 
where F (z) and G (z) are given by (17) 



h 



G (z)j 'V3 [»!(«) 



■^i(z) 



(31) 



6. A More Tractable Solution to (1) 

The general solution to the singular integral equation (8) is now readily obtained from (31). 
By the Plemelj formulae (13b), the general solution to (8) is given by 



4>{x)=F + {x)-F_{x) 



h 



=a (Fo + (x)-Fo-(x))+^r-(F u (x)-F u (x)) 
*V3 



(32) 
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A_(*^)],i 



=a (t>o(x)-{ — 7^1-1=- j'lOO 

where <t> (x) is given by (15) and 4> (x) is given by (24). 

In turn, by the relations (9), (15), and (22a), the general solution \p(x) to (1) is given by 

hoc ^ 
Hx)=a TJTk (x)+-^ M*), 0<s<l, (33) 

where k (x) is given by (3a) and \p (x) is given by (23). This expression for \[/(x) is certainly a 
more tractable one than the expression for \(/{x) given by (3). 

This result can now be used to finish the analysis of the wave-guide problem. An identi- 
fication in (2c) says that 

E'(y)=\P(cos (y))=a ^Tk (cos (y))H = cos (?y)^ (cos (y)). (34) 

Let 

Vtt& (cos (?/))=^o(cos (?/)) 
so that 

£"(7/)=a iAo(cos (y))-\ = cos (y)M<x>s (y)). (34a) 

7rV3 

Following the analysis of Bierman [1], E'(y) is odd and continuous at y=0. Therefore, 
E"(0) = 0, or equivalently 

lim ( a ^o(cos (y))+—= cos (y)$ (coa(y)))=0. (35) 

v^o \ 7rV3 / 

Relation (35) thus determines a uniquely, 

= -h 
tV3 

To obtain E(y), eq (34) is integrated, and the constant of integration is chosen so that 

as required in [1]. 

E(y)=-^j= f'* (cos (a)^o(cos («r))-*,(00B (a)))«fcr. (37) 

7rV3 Jy 

The reflection coefficient, p, is obtained from (2a). 



2ttV3 



-j=ikE (l-\-3p) (cos (o-)^o(cos(o-))— ^ (cos(o-)))(T(Za-. 
/3 Jo 
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Hence 



where 



1+Af 
p 1-3M 



ik C* 12 ( m \ 

M= p Vcos (o-)^o(cos (a))— \l/ (cos (cr))J<rdo 

27tV3Jo 



rV3. 
Once p is computed, T can be computed by (2b). 
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